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Abstract. We explore the non-equilibrium escape dynamics of long-range interacting 
ions in one-dimensional traps. The phase space of the few ion setup and its impact 
on the escape properties are studied. As a main result we show that an instantaneous 
reduction of the trap's potential depth leads to the synchronized emission of a sequence 
of ion pairs if the initial configurations are close to the crystalline ionic configuration. 
The corresponding time- intervals of the consecutive pair emission as well as the number 
of emitted pairs can be tuned by changing the final trap depth. Correlations between 
the escape times and kinetic energies of the ions are observed and analyzed. 
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1. Introduction 

The controlled manipulation of ions in traps has triggered several areas in physics. 
Especially the impact on the development of high precision measurements in 
spectroscopy [1] and on concepts for implementing quantum simulations (see [21 [H H] 
and references therein) has been substantial. An interesting effect observed in ion- 
trap resonators is the synchronization of ionic motion [5] with applications in mass 
spectrometry. Furthermore, the study of the stability and the loss processes of ions 
shows that in linear electrostatic traps up to 30 % of trapping efficiency is reachable 
[Hj. By applying the (5-kick method a compression of the velocity spread of stored ions 
in these traps has been demonstrated [7], which can be exploited for cooling. Another 
widely-used ion trap is the so-called Paul trap [8] which has been used to observe and 
monitor single ions [HI HDl E] or to study the stability of large Coulomb ion crystals 
[T2| [T3| [H] . Originally, the first experimental observation of a phase transition from a 
cloud of ions to a crystalline structure has been reported in Ref. [I5]. In this work, the 
ions were cooled by means of a light background gas and by varying the trap potential 
and the cooling force, repeated crystallizations and meltings could be monitored. Phase 
transitions have been observed for laser-cooled Mg"*" as well as Hg+ ions in a Paul trap 
in Refs. [16], and [17], respectively. The phase transition of a two ion crystal to a 
cloud has been explained by the nonlinear aspects of the Paul trap, i.e. the melting 
of a crystal has been interpreted as an order to chaos transition induced by ion-ion 
collisions [T8l [T9l [201 [2T] . Further research on nonlinear phenomena in trapped ion 
systems dealt with transient chaotic regimes preceding the transition point of the Paul 
trap [22], frequency-locked motion of two ions in the Paul trap [23] and with nonlinear 
effects in the dynamic Kingdon trap [21]. In Ref. [251(261 [27] the dynamics of Coulomb 
chains with many ions in a harmonic potential has been studied. The authors evaluate 
analytically the long- wavelength modes and the density of states in the short- wavelength 
limit. Using these results they are able to derive the critical transverse frequency of 
the trap below which the chain is unstable with respect to excitations of transverse 
vibrations. 

In this work we study the escape dynamics of a chain of ions in a finite box-like 
smooth potential. The non-equilibrium properties of this one-dimensional system of 
long-range interacting ions will be explored for the case when the trapping potential is 
instantaneously reduced (we will denote this in the following by a "ramp" of the potential 
well) such that the emission of particles from the trap is possible. Our aim is to unravel 
and analyze the non-equilibrium dynamics as well as the escape properties of ions from 
such a one-dimensional potential well. We find that above an energetic threshold the 
energy shell reaches out to infinity thus making the escape of ions from the potential well 
possible. The emission process is dominated by the unstable periodic orbits and their 
manifolds in the sense that they lead to dominant peaks in the distribution of the kinetic 
energy of the escaping ions. Ramping up the potential we observe a strong correlation 
between the escape time and the kinetic energy of the emitted ions. Most importantly 
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we show that the ionic emission occurs in a sequence of synchronized pairs whose 
properties are tunable via the trap parameters. In particular by a ramp of the well to a 
corresponding depth the number of emitted pairs as well as the time difference between 
consecutively emitted pairs and their kinetic energies can be adjusted. Accordingly, this 
escape process from the trap yields a deterministic source of correlated ion pairs with 
selectable dynamical properties. Such a source could be used for performing scattering 
experiments in which one needs energetically and spatially correlated particles. 

This paper is organized as follows: In Sec. |2] we introduce the Hamiltonian and 
describe the corresponding scaling transformations. Sec. [3] is devoted to the study of 
the classical phase space of a single and two ions in the static potential well. In Sec. H] 
the escape dynamics of many ions from the trap after a ramp of the potential depth is 
discussed. Finally, our summary is provided in Sec. [5l 



^(^) = -^^^^ (2) 



2. Setup 

The Hamiltonian (in Gaussian units) describing the dynamics of n Coulomb-interacting 
particles with equal charge q and mass m in an external potential V{x) is given by 

n 2 " 1 2 

^-^ 2m ^-^ 2 ^-^ ^-^ \Xi — Xi 

i=l 4=1 i=l jyti ' ■> ' 

where x = (xi, . . . , p = {pi, . . . ,pn) are the positions and the momenta of the 
ions, respectively. The ions are numbered according to their alignment, i.e. xi is the 
position of the leftmost ion and pi its momentum etc. For V{x) we choose the following 
functional form 

1 + {cxf 

which corresponds to a smooth potential well of depth Vq. c fixes the width of the 
well via V{l/c) = — Vo/2. Fig. [1] (a) shows the potential V{x). For this functional 
form of the potential the Coulomb crystal configuration is an almost equidistantly 
spaced ion chain. The one-dimensional confinement of the dynamics can be achieved for 
example in RF-traps by appropriately adjusting the trap frequencies along the transverse 
dimensions, i.e. when the kinetic energy of all ions is small compared to the transversal 
harmonic confinement frequency [251 ISSl [27] . Hence, the dynamics of the system is one- 
dimensional. Apparently, the parameter space of the Hamiltonian is four-dimensional: 
the depth of the well Vq, the width of the well determined by ^, the mass of the ions m 
and the charge of the ions q. It is therefore useful to apply a scaling transformation to 
dimensionless variables: 

X ^ X = ex, p ^ p = — p 

m 

2 

If the Hamiltonian is scaled additionally according to H ^ H' = the result reads 

H'{x'j) = yp^-y^^ + -yy^^^, (3) 

1 = 1 * J=l Jytl ' ' J ' 
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Figure 1. (a) Potential V{x) as given by Eq. 
bounded dynamics of a single ion in V{x). 



([2]). (b) Phase space portrait for 



with a scaled charge q' = \J and potential Vq = ^Vq. Obviously, this transformation 
yields a reduction of the dimension of the parameter space to two: the potential depth 
Vq and the charge q' . For reasons of clarity the prime is omitted in the following. The 
corresponding Hamiltonian equations of motion 
dH . OH 

Xi = Pi = (4) 

opi dxi 

are solved numerically using the Burlisch-Stoer algorithm [28]. 



3. Phase space for few ion systems 

Let us now turn to the discussion of the phase spaces for setups with few ions in the 
potential well. All results in this section are obtained for q = 1 and Vq = 5. Since 
the energy is a constant of the motion the one ion system is integrable and the two 
dimensional phase space can be visualized straightforwardly. Fig. [1] (b) shows the level 
curves for bounded oscillatory motion which are simply the contour lines of the energy. 
Obviously, the initial momenta for these trajectories must obey p < \/2Vq. For two ions 
in the trap a representative visualization is still possible. Due to energy conservation 
the dynamics in the four dimensional phase space is confined to a three-dimensional 
energy shell. The lowest energy state of two ions in the static well is the so-called 
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Coulomb crystal and corresponds to a single point in phase space. For our choice of the 
parameters for Vq and q this configuration possesses the energy Ec = —9.0126. The 
energy shells in phase space with E > Ec are visualized by the corresponding Poincare 
surfaces of section (PSS) using the condition pi = 0. The PSS is made such that we 
report crossings from pi > to pi < 0, i.e. we plot the momentum against the position 
of the right (second) ion, when the left one reaches its inner turning point. Fig. |2] shows 
the PSS for different values of the total energy. For E = —9 the dynamics is completely 
integrable (Fig. |2] (a)). At Si = (xs = 0.537,^2 = 0) and ^2 = (xs = 0.606, p2 = 0) 
there are stable period one orbits, which correspond to the two normal modes, i.e. in 
the case of 5*1 both ions oscillate with opposite phase such that the center of mass is 
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at rest, whereas 5*2 belongs to the center- of -mass mode (oscillation in phase). With 
increasing energy the dynamics becomes richer. Elliptic islands embedded in a chaotic 
sea appear (Fig. [2] (b)), whose area in the PSS grows until E = —5. For even larger 
energies one observes a growth of the central integrable part. For E = —4.59 initial 
conditions in the chaotic able to leave the well, i.e. the energy shell reaches in 

the Xi, a;2-direction to infinity which can be seen as an excrescence on the right hand 
side in Fig. [2] (c). For higher initial energies the size of the opening of the phase space 
towards infinity gets bigger which allows for a quicker escape of the particles from the 
chaotic sea to infinity, e.g. in Fig. [2](d) the initial conditions outside of the shown area of 
closed curves belong to trajectories where the ions leave the trap rapidly. Furthermore, 
with increasing energy the remaining integrable part shrinks, until for E > —3.6 no 
integrable structure in phase space is present. In the energy range E > —3.6 unstable 
periodic orbits and their manifolds organize the dynamics in phase space. In order to 
detect them we apply a two-dimensional Newton-Raphson-method [28] in the PSS. Fig. 
|3] shows the stable and unstable manifolds of an unstable period one orbit located at 
U = (x2 = 0.063, p2 = 0) for E = —2. The fiow of the unstable manifolds is constructed 
by distributing a whole line segment of initial conditions along the directions of the 
eigenvectors of the underlying monodromy matrix. Both manifolds reach in the x- 
direction out to infinity. In Fig. [3] this is only partially visible because the part of the 
manifolds which reach out to infinity is cut off in the x-range of the well. Below we will 
see that the manifolds are important for the understanding of the escape of ions from 
the well. 

3.1. Escape of a single ion from the well 

In this section the properties of the single ion escape process from the well are discussed, 
i.e. one ion possesses via the Coulomb interaction with its partner sufficient energy 
such that it can be emitted. We do not consider the two ion escape processes which 
is possible for E > 0.42 because in this high-energy regime the immediate escape of 
both ions dominates and no new time scales or processes emerge. In order to avoid 
artificial transient effects we have to ensure that the initial conditions are distributed 
uniformly over the energy shell. To this end, we briefly describe in the following 
the preparation of the initial ensemble. First, the positions of the ions (xi,X2) are 
distributed randomly in the interior of the well, which we deflne to be \x\ < 2.5. 
Due to energy conservation, the corresponding momenta lie on a circle with radius 
r{xi,X2) = \/p\+ pI = \/lm{E — 'V{x\,X2)) for each pair [xi.x-i). On these circles in 
momentum space the points are distributed with the same density as in position space. 
Altogether, this procedure yields a uniform distribution of quadruples {x\,X2,P\,P2) on 
a given energy shell. From the previous discussion of the phase space properties it is 
known that unbounded motion is only possible for energies E > —4.59. Accordingly, 
we distribute approximately 10^ initial conditions uniformly over different energy shells 
with E > —4.59 and record the escape time and the kinetic energy at the border of 




the well (note that all observed effects and phenomena are qualitatively not affected by 
the position of the detector, where the escape time and kinetic energy are recorded). 
Generally, with increasing initial energy the region in phase space allowing for escape 
grows and the particles leave the well more quickly. Nevertheless, we observe also 
events where it takes long time intervals until emission occurs. Fig. H] shows the escape 
times on a logarithmic scale for initial conditions on the energy shell E = —2 which 
lie additionally in the plane pi = formed by the PSS. The border of the grey area 
defines the border of the energy shell. A comparison of Fig. |4] with Fig. [3] (a) shows 
that the longest escape times belong to initial conditions close to the stable manifold 
of the period one orbit U. This orbit is marginally unstable because it originates from 
Si which is one of the dominant stable periodic orbits. Consequently, the particles 
starting in the vicinity of the stable manifold of U approach the fixed point, stay in 
its vicinity for a long time and escape afterward via the unstable manifold leading to 
infinity. Fig. [5] (a) shows the probability density of the kinetic energy of the escaping 
ion p{E]^i^) for initial energy E = —4, where ii^kin is recorded at the border of the well 
|x| = 2.5. p(-E'kin) possesses a sharp peak at the maximal value for the kinetic energy, 
which is determined by energy conservation. Obviously, the Coulomb term as well as 
the potential energy of the escaped ion is zero. Furthermore, the minimum energy of 
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Figure 4. Logarithm of the escape time logt for initial conditions in the PSS pi = 
on the energy sheh E = —2. 

the residual ion is = —Vq. Thus the maximum kinetic energy of the escaped ion 
is -Ekin.max = E + Vq which amouuts in our case to -Ekin.max = 1- Since there is still a 
remaining Coulomb contribution (the detector is at a finite distance from the trap) for 
the distribution shown in Fig. [5] (a), the maximum kinetic energy -Ekin.max is smaller, i.e. 
p(£^kin) has a sharp cut-off at -Ekin = 0.65 and not at .Ekin = 1-0. The general appearance 
of p(i?kin); namely an increase followed by a plateau-like behaviour, can be understood 
intuitively. For E = —4 the phase space contains a regular island of bounded dynamics 
in the potential well (see Fig. [2] (d)). Consequently, the ions starting in this integrable 
structure do not contribute to the escape. Due to the fact that these initial conditions 
correspond to small kinetic energies, the distribution p(i?kin) decreases for £^kin ~^ 0. Fig. 
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Figure 5. Probability distribution of kinetic energies (a) and escape times (b) of the 
escaping ion for the energy shell E = —4 and Vb = 5. 



M (b) shows the distribution of escape times p{t) which has a dominant peak followed 
by a steep drop. With decreasing initial energy the tail of the distribution becomes 
more pronounced until a collapse occurs for E = —4.6 and no escape of ions is possible 
anymore. The peak of p{t) corresponds to the immediate emission process. A further 
analysis of the correlation between the kinetic energy (Fig. [5] (a)) and the escape time 
(Fig. |5](b)) reveals that the peaks in p(£'kin) belong to initial conditions with particularly 
long escape times and thus to particles starting in the vicinity of the stable manifolds 
of marginally unstable fixed points (see Fig. H] and the corresponding discussion of the 
marginally unstable periodic orbits). The width of these peaks of p{Ekin) corresponds 
to the volume in phase space which gets transported along the unstable manifold out of 
the trap. The structure of these manifolds becomes more complicated with decreasing 
energy leading to a sequence of peaks in p(i?kin). Finally, the channel of the energy 
shell towards infinity becomes thinner. Accordingly, the distribution of kinetic energies 
develops into a delta distribution located at i^^in = 0.4 for E — —4.6. 



4. Escape dynamics after a ramp of the well 

In this central section of our work the dynamics following a ramp of the well with many 
ions is explored. A ramp means that the depth of the well is reduced instantaneously. 
Due to the fact that the number of degrees of freedom increases according to the number 
of particles, the energy shells of the static system for many ions are topologically 
complicated objects of high dimensionality. Hence a uniform distribution of initial 
conditions on an energy shell is not only difficult to achieve but in particular also not 
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a natural starting point for an experimental setup. Instead we use as initial conditions 
an ensemble which is sampled in position space randomly in small intervals Ax around 
the Coulomb crystal configuration with all initial momenta set to zero. For the interval 
widths Ax we choose, as a typical value, the twentieth part of the mean crystal distance. 
Instead of the uniform distribution on a given energy shell this is indeed a more realistic 
scheme with respect to an experimental preparation of such a setup. For the charge of 
the ions q and the depth Vq we choose g = 1 and Vq = 30 such that the Coulomb crystal 
exists for up to ten ions in the well. In a second step the well is ramped up by AV^ and 
the resulting dynamics is observed and analyzed. 

Three ions For AV < 20 the ions oscillate with comparatively small amplitudes and 
no significant energy transfer among them occurs. However, for a further reduction of 
the well's depth Vq the Coulomb interaction gains importance such that the dynamics 
becomes more and more chaotic: for a larger AV^ the ions oscillate with larger 
amplitudes and thus the relative distance between them can become small. Therefore, 
the probability increases that in the course of the dynamics after the ramp an ion gains 
enough kinetic energy to be emitted. Still, due to energy conservation there is an energy 
threshold for a single ion-escape to become possible. Nevertheless, above this threshold 
not every initial condition leads to an emission-event. For Vq = 2.3 the escape of a 
single ion becomes energetically possible. For this energy which is slightly above the 
mentioned threshold value the escape process is typically characterized by a phase of 
motion during which two neighbouring ions move in the same direction. This leads to a 
strong enhancement of the kinetic energy of the third outermost ion which gets "kicked" 
out of the potential well. After the escape of one ion the relevance of the Coulomb 
repulsion decreases and the remaining ions start quasi-periodic oscillations around the 
equilibrium positions of the two-ion Coulomb crystal. In Fig. [6] the distribution of the 
escape times for a ramp to Vq = 2.1 is presented, which shows asymptotically (large t) 
an algebraic decay. The long tail is due to initial conditions with total energies hardly 
above the energy threshold. Consequently, for these events the energy exchange process 
leading to an ion-escape has to be very efficient meaning that two ions must transfer a 
large fraction of their kinetic energy to the third ion such that the latter can leave the 
well. Accordingly, a larger AV of the ramp of the potential does not necessarily yield a 
higher kinetic energy of an individual escaping ion because the emission occurs also in 
case a smaller portion of the kinetic energy of two ions is transfered to the third one. By 
a ramp of the well to even smaller depths Vq < 1.5 the two ion emission process becomes 
possible. Finally, for Vq < 1.3 all initial conditions yield two ion escape. In this regime 
the emission occurs generally in a synchronized manner, i.e. as pairs where the difference 
between the escape times Atp = ti — ^2 of the two ions of a pair is small compared to all 
other time scales (see below for a more detailed discussion of this process). We remark 
that by fine-tuning AV it is also possible to achieve the situation where Atp can become 
large. Such a two-ion emission occurs if the ramping of the well AV is tuned to be only 
slightly above the energy threshold for two-ion emission. These events are rare and their 
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Figure 6. Probability distribution of escape times p{t) of a single ion when three ions 
arc loaded into the trap and the well is ramped from Vq = 30 to Vq = 2.1. 

fraction remains below 0.75 % in the regime 1.5 > Vq > 1.3. 

Ten ions For ten ions the Coulomb interaction is stronger and thus we expect a larger 
impact of the instantaneous reduction of the potential depth. Fig. [7] shows the dynamics 
of trajectories for a ramp of the trap from Vq = 30 to Vq = 2. The single initial condition 
used for the simulation corresponds to the configuration of the Coulomb crystal for 
Vq = 30. In this case we observe that eight ions leave the well. Obviously, the emission 
of the ions happens in a sequence of four pairs. Every time a pair is emitted the 
corresponding two ions leave the trap in opposite directions. 

In the following we study the robustness of this behaviour with respect to the 
choice of the initial ensemble, i.e. for a random sampling in the vicinity of the Coulomb 
crystal. Accordingly, we take 10^ different initial ensembles of ions, ramp up the trap 
by AV and record the kinetic energy as well as the escape time of the emitted ions. 
Thus all following results are ensemble properties. Fig. [8] shows the logarithm of the 
frequency distribution log(p(£'kin, i)) of the kinetic energy and the escape time for a 
ramp from Vq = 30 to Vq = 4. In this case six ions escape from the trap. Each 
of the three bright dots in Fig. |8] corresponds to two emitted ions. In the inset the 
frequency distribution of the time differences p(Atp) between the escape of the first and 
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Figure 7. Trajectories of the ten ions when the potential depth is instantaneously 
reduced to Vq = 2. For Vq = 30 the initial condition corresponds to the Coulomb- 
crystal. 

the second ion and the fifth and the sixth ion is shown. p{Atp) belonging to the third 
and fourth ion is almost identical to the corresponding distribution for the first and 
the second ion. Obviously, Atp is very small compared to the interpair emission time 
intervals Atpp. Additionally, the distribution p(Atp) decreases rapidly, i.e. to a very good 
approximation the emission process can indeed be regarded as pairwise. An intuitive 
physical explanation of this result is given in the following: Before the ramp the initial 
ensembles are spatially close to the highly symmetric Coulomb crystal configuration. 
Due to the instantaneous reduction of the potential depth much potential energy is 
transfered into the system. Subsequently, these now excited chains relax by the emission 
of ions. Since the symmetry of the trap is preserved by the ramp, the forces acting on 
the outermost ions immediately after the instantaneous reduction of the potential have 
the same absolute value but opposite signs and consequently the ions escape in opposite 
directions within a small time-interval. Furthermore, Fig. |H] reveals that Atpp increases 
within a sequence of pairs whereas the kinetic energy decreases. This behaviour can be 
explained as follows. The force which pushes the ions outwards is largest for the first 
pair and decreases with every subsequently emitted pair. Additionally, every escaping 
pair is slowed down by the Coulomb repulsion of the previously emitted ions. In Fig. 
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Figure 8. (Color online) Logarithm of the frequency distribution of the kinetic energy 
and the escape time log(p(i5kin, 0) foi" ^ ramp from Vq = 30 to Vq = 4, where six 
ions leave the well. Each white spot belongs to two emitted ions. In the inset the 
distribution of the differences of the escape times p{Atp) between the first and the 
second ion (solid line) and the fifth and the sixth ion (dashed line) are shown. 



[7] these two characteristics can be seen as a "fanning out" of the trajectories with a 
decreasing slope. Thus the escape velocity becomes smaller with every emitted pair. By 
ramping the trap to different potential depths in the regime where the emission of six 
ions occurs, the time difference between the consecutive pairs Atpp can be tuned, i.e. 
Atpp decreases with increasing AV. Moreover, since the total energy of the ensemble 
after the ramp is higher for larger AV, the kinetic energy of each emitted pair increases 
as well. 

So far we have discussed the situation when the ramp leads to an emission of an 
even number of ions. Yet, it is also possible that after the ramp only an odd number 
of ions can leave the trap. For these cases the scenario is the following. First one 
observes the emission of a maximal even number of ions in terms of synchronized pairs 
according to the discussion provided above. Afterwards, a single ion leaves the well. 
Its escape from the trap can be delayed significantly with respect to the emission of 
the last pair. Such an event occurs when the ramp is chosen only slightly above the 
necessary energy threshold and can be explained straightforwardly. After the emission 
of the pairs the remaining ions oscillate chaotically in the trap. Accordingly, it can 
take a very long time until, in the course of the dynamics, the randomly distributed 
kinetic energy accumulates onto a single ion which can then be emitted. Consequently, 
the distribution of the kinetic energies of the last single escaping ion is not peaked at a 
certain value as for the pairwise emission but it shows an exponentially decreasing tail. 
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Figure 9. (Color online) Logarithm of the frequency distribution of the kinetic energy 
and the escape time log(p(iJkin, t)) for a ramp from Vq = 30 to Vq = 4 of an asymmetric 
well. In this case four ions are emitted from the trap. The inset shows the potential 
described by Eq. ([5|) for Vo = 30 and ck = ^ 



Let US now discuss the symmetry robustness of the synchronized pair emission 
process. To this end we modify the functional form of the potential by multiplying it 
with an asymmetric exponential term, i.e. 



In Fig. |9]the modified potential Eq. ([5]) is shown in the inset for Vq = 30 and « = ^• 
For X — )■ — oo we have V{x) — t- — oo. However on the scale on which the dynamics is 
observed this behaviour is suppressed by the large power of the denominator and the 
fact of a small value of the pref actor a. Due to the parity asymmetry of the potential 
the Coulomb crystal and thus the sampling of the initial ensemble is also asymmetric 
with respect to x = 0. Fig. [9] shows log(p(i?kin, i)) for a ramp from Vq = 30 to Vq = 9. 
In this case four ions escape from the trap. The first ion leaves the well in positive 
x-direction at t ^ 0.52 with kinetic energy E]^[^ ^ 9.59 (see Fig. [9]). The second ion 
which has less kinetic energy E'kin ~ 9.22 escapes to — oo at t ~ 0.56. Afterwards, the 
emission of a pair with E'kin ~ 3.68 occurs at t ^ 1.02. This behaviour can be traced 
back to the shape of the trap. Obviously, due to the asymmetry of the potential the 
leftmost ion has a lower potential energy. Thus compared to the rightmost ion a larger 
fraction of its kinetic energy is needed for the emission from the trap after the ramp of 
the potential depth. Since the potential difference due to the asymmetry of the trap 
decreases towards the origin, the difference of potential energies is less for the next two 
outer ions. Thus they escape pairwise again. For a larger ramping AV one observes also 




(5) 
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an energetic and temporal separation of the emission for the later escape events, that is 
for an appropriate choice of the third and fourth ion etc can be separated as well. 
Due to the fact that these results are obtained for a random sampling of the ions' initial 
position in the vicinity of the corresponding Coulomb crystal we expect that effects 
like micromotion do not alter qualitatively the results. Furthermore, this statement is 
supported by our numerical observation that no significant changes occur when the ions 
are additionally given a small initial momentum, i.e. the emission properties of the ions 
from the trap are very robust. 

5. Summary 

We have studied the nonlinear escape dynamics of long-range Coulomb interacting ions 
from a one-dimensional trap. As a necessary prerequisite we have first accounted for the 
few ion setup, its phase space structure, and observed that the phase space of two ions 
in the static potential well can be regular, mixed or chaotic depending on the energy. 
For energies larger than a certain threshold the energy shell extends towards infinity 
and consequently the escape of ions from the trap is possible. It has been elucidated 
that the manifolds of unstable periodic orbits play a crucial role in this emission process. 
Particles starting in the vicinity of the stable manifold show the longest escape times 
from the well and contribute to dominant peaks in the distribution of the kinetic energy 
of the escaping ions. 

Our main result concerns the escape dynamics of an ensemble of ions initially in a 
configuration close to the Coulomb crystal when the potential depth is decreased by an 
instantaneous ramp of the trap. We have observed that the emission of an even number 
of ions occurs in a sequence of pairs. Accordingly, every time a pair is emitted the 
corresponding two ions escape in opposite directions. A deeper analysis of the escape 
time and the kinetic energies reveals a strong correlation between both quantities, that 
is each pair leaves the trap with a definite kinetic energy. By adjusting the value 
of the ramp the number of the emitted pairs, the time differences between the 
consecutive pair emission events and the average kinetic energy of each pair can be 
varied. If the ramp is chosen such that an odd number of ions can escape, one observes 
after a maximal even number of ions escaping pairwise the emission of a single ion. Its 
escape can be delayed significantly. Moreover, the kinetic energy distribution possesses 
an exponentially decreasing tail, because the dynamics leading to its escape is chaotic. 
Furthermore, we have studied the symmetry robustness of the escape properties. It has 
been found that parity asymmetry of the potential leads to a temporal and energetic 
separation of the pairwise emission. We have also considered the robustness of our 
results with respect to variations of the potential by choosing a less steep and box-like 
potential V[x) and found that this does not alter qualitatively the results. 

Obviously, the ramp of the trap from the regime where the initial ensemble is 
confined in a configuration close to the Coulomb crystal to a lower potential depth 
can be exploited to create a temporally synchronized source of ions. By appropriately 
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choosing AV a desired number of emitted ions (pairs) and their escape properties hke 
escape time and kinetic energy can be adjusted. 

Acknowledgments 

The authors acknowledge gratefully helpful discussions with M. Drewsen and F. K. 
Diakonos. This work has been performed within the Excellence Cluster Frontiers in 
Quantum Photon science, which is supported by the Joachim Herz Stiftung. Financial 
support by the DAAD in the framework of an exchange program with Greece is 
acknowledged. Part of this work was financially supported in the framework of 
the excellence initiative of the federal and state government of Germany, i.e., the 
FRONTIER program of the university of Heidelberg. 

References 

[I] R. C. Thompson, Adv. At. Mol. Phys. 31, 63 (1993). 

[2] D. Leibfried, R. Blatt, C. Monroe, and D. Wineland, Rev. Mod. Phys. 75, 281 (2003). 
[3] H. Haffner, C. F. Roos, and R. Blatt, Phys. Rep. 469, 155 (2008). 

[4] M. Johanning, A. Braun, N. Timoney, V. Elnian, W. Neuhauser, and C. Wunderhch, Phys. Rev. 
Lett. 102, 073004 (2009). 

[5] H. B. Pedersen, D. Strasser, S. Ring, O. Heber, M. L. Rappaport, Y. Rudich, I. Sagi, and D. 

Zajfman, Phys. Rev. Lett. 87, 055001 (2001). 
[6] H. B. Pedersen, D. Strasser, O. Heber, M. L. Rappaport, and D. Zajfman, Phys. Rev. A 65, 042703 

(2002). 

[7] S. Goldberg, D. Strasser, O. Heber, M. L. Rappaport, A. Diner, and D. Zajfman, Phys. Rev. A 68, 

043410 (2003). 
[8] W. Paul, Rev. Mod. Phys. 62, 531 (1990). 

[9] W. Neuhauser, M. Hohenstatt, P. Toschck, and H. Dehmelt, Phys. Rev. Lett. 41, 233 (1978). 
[10] W. Neuhauser, M. Hohenstatt, P. Toschek, and H. Dehmelt, Appl. Phys. 17, 123 (1978). 

[II] W. Neuhauser, M. Hohenstatt, P. Toschek, and H. Dehmch, Phys. Rev. A 22, 1137 (1980). 

[12] M. Drewsen, C. Brodersen, L. Horneka^r, J. S. Hangst, and J. P. Schifffcr, Phys. Rev. Lett. 81, 
2878 (1998). 

[13] N. Kjffirgaard, K. M0lhave, and M. Drewsen, Phys. Rev. E 66, 015401 (2002). 
[14] N. Kja;rgaard, and M. Drewsen, Phys. Rev. Lett. 91, 095002 (2003). 
[15] R. F. Wuerker, H. Shelton, and R. V. Langmuir, J. Appl. Phys. 30, 342 (1959). 
[16] F. Diedrich, E. Peik, J. M. Chen, W. Quint, and H. Walther, Phys. Rev. Lett. 59, 2931 (1987). 
[17] D. J. Wineland, J. C. Bergquist, W. M. Itano, J. J. Bollinger, and C. H. Manncy, Phys. Rev. Lett. 
59, 2935 (1987). 

[18] R. Bliimel, J. M. Chen, E. Peik, W. Quint, W. Schleich, and H. WaUher, Nature 334, 309 (1988). 

[19] J. Hoffnagle, R. G. DcVoc, L. Rcyna, and R. G. Brewer, Phys. Rev. Lett. 61, 255 (1988). 

[20] R. Bliimel, C. Kapplcr, W. Quint, and W. Wahher, Phys. Rev. A 40, 808 (1989). 

[21] R. G. Brewer, J. Hoffnagle, R. G. DeVoe, L. Reyna, and W. Henshaw, Nature 344, 305 (1990). 

[22] R. G. Brewer, J. Hoffnagle, and R. G. DeVoe, Phys. Rev. Lett. 65, 2619 (1990). 

[23] J. Hoffnagle, and R. G. Brewer, Phys. Rev. Lett. 71, 1828 (1993). 

[24] R. Bliimel, Physica Scripta T 59, 369 (1995). 

[25] G. Morigi, and S. Fishman, Phys. Rev. Lett. 93, 170602 (2004). 

[26] G. Morigi, and S. Fishman, Phys. Rev. E 70, 066141 (2004). 

[27] G. Morigi, and S. Fishman, J. Phys. B: At. Mol. Opt. Phys. 39, 221 (2006). 



Correlations and pair emission in the escape dynamics of ions from one- dimensional trapsll 
[28] W. H. Press, Numerical recipes in C, Cambridge University Press, (1992). 



